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An analytical model has been developed to simulate the material dif- 
ferentiation of a cometary nucleus composed of water ice, putative 
unclathrated CO^ ice and silicate dust in specified proportions. Selective 
sublimation of any free CO. ice present in a new comet would produce a sur- 
face layer of water ice and dust overlying the original CO rich material. 
This surface layer reduces the temperature of buried CO^ ice and restricts 
the outflow of gaseous CO^,. On each orbit, water sublimation at smaller 
heliocentric distances temporarily reduces the thickness of the water ice 
and dust layer and liberates dust. Most of the dust is blown off the 
nucleus, but a small amount of residual dust remains on the surface. (cf. 
Houpis et al., 1985) Our model includes the effects of nucleus rotation, 
arbitrary orientation of the rotation axis, latitude, heat conduction into 
the deep interior of the nucleus and restriction of 00 ^ gas outflow by the 
water ice and dust layer, features that were not included in the Houpis et 
al. model. Specifically, we investigate the effects of the permeability of 
the surface water ice layer, the nucleus rotation rate and the latitude. 
The loss rate of CQ„ and the resultant depth of CO^ ice are shown to be 
most sensitive to the permeability of the water ice and dust layer. For a 

homogeneous, initially unmantled comet placed in the orbit of comet Halley, 

it is shown that the CO- ice attains a steady state or cyclic relationship 
between CO^ depth and orbital position within several revolutions. If the 
nucleus contains 25% by mass of (XL, our results indicate that CO^ ice is 
always within several meters from the surface at axrp location for a nucleus 

of low obliquity and that CO^ iee is nearest to the surface at the equator 

shortly after perihelion. Under these conditions the sublimation of CO- 
ice is always significant and becomes the dominant gaseous species beyond 
4AU, This result is probably generally valid for unmantled portions of 
most comets and qualitatively simulates the behavior of an abundant, highly 
volatile component in an H 2 O/silicate matrix. Comparison of these and simi- 
lar results with observations could yield information regarding the permea- 
bility and chemical composition of cometary material and suggest sampling 
strategies to minimize fractionation effects. 
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